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While the diagonalization of a quadratic bosonic form can always be 
done using a Bogoliubov transformation, the practical implementation for 
systems with a large number of different bosons is a tedious analytical task. 
Here we use the coupled cluster method (CCM) to exactly diagonalise such 
complicated quadratic forms. This yields to a straightforward algorithm 
which can easily be implemented using computer algebra even for a large 
number of different bosons. We apply this method on a Heisenberg system 
with two interpenetrating square lattice antiferromagnets, which is a model 
for the quasi 2D antiferromagnet Ba2Cu304Cl2- Using a four-magnon spin 
wave approximation we get a complicated Hamiltonian with four different 
bosons, which is treated with CCM. Results are presented for magnetic 
ground state correlations. 

PACS numbers: 75.50.Ee,75.30.Ds,31.15.Dv,42.50.Ls 



1. Introduction The Model 

It is always possible to diagonalize quadratic bosonic forms (which appear 
frequently in physics) using a Bogoliubov transformation ||, but it is a tedious 
analytical task to find one for a complicated form with many different magnons. 
Therefore we want to show here how the coupled cluster method (CCM), one of 
the most powerful and universal techniques in quantum many-body theory (s. (^] 
and references therein), can be used in a straightforward scheme to find the exact 
ground state of such a form. 

To be concrete we consider the 2D spin 1/2 Heisenberg model 

H = Jaa ^ SjSj + Jbb ^ SjSj + Jab ^ SjSj, (1) 

{i£A 1 ,j£A 2 ) {ieB 1:] eB 2 ) (ieAjGB) 



(1) 



2 



which is related to the situation in Ba2Cu304Cl2 ||, f|, a layered quantum an- 
tifcrromagnet showing significant differences to its parent cuprats like La2CuC>4 
(see e.g. jl]] for recent experiments). In contrast to La2Cu04 we have two different 
types of Cu-sites in the Cu-O-planes, namely there are additional Cu(B) atoms 
located in the centre of every second Cu(A)-0 2 square. Within the Cu(A) sub- 
system we have a strong 180° Cu-O-Cu superexchange yielding to strong antifer- 
romagnetic couplings (Jaa) between Cu(A) atoms, whereas the couplings within 
the Cu(B) subsystem (Jbb) and between the subsystems (Jab) are weaker. A 
recent calculation of Jaa, Jbb f§], finding Jaa ~ IOJbb (both antiferromagnetic) 
agrees with the experimental values Q. There are also some arguments Q for a 
ferromagnetic \Jab \ ~ Jbb- 

In the classical ground state (0) shows for \Jab\ < 2\J JaaJbb a Neel like 
order for the two subsystems A and B, where the energy is degenerated with 
respect to the angle ip between the spins of these two subsystems. 

2. The Method 

In this paper we study the ground state properties of (|l|), using a four- 
magnon linear spin wave approximation |j) around the classical ground state, 
i.e. for each of the four sublattices A\ , A 2 , B\ , B 2 of the two coupled bipartite 
antiferromagnetic square lattices we introduce different bosonic operators. Thus 
we get for (|l|) 

2N 

H = — -s 2 (2J A a + Jbb)+Y, H ^ with ( 2 ) 
6 k 

#k = 4 Jaas (af k aik + a^ k a 2 k - JkAA [a^aj-k + a ik«2-k] ) 

+ 2J B BS (fr/" k &lk + &2k^2k - 7kSB [K^2-k + °lk02-k] ) 

+Jabs(1 + cosp)/2 (foj^aik + h^a^y. - &2k a ^-k _ &2kGi-k) ll A B 

+ JABS (1 - COS^) /2 (&2k a 2k + &2ka2k _ b rk a 2-k ~ °lka2-k) 7kAB 

+Jabs (1 - cos^) /2 (&2k a ik + hka^ - &i k a^_ k - 6ikai-k) 7kAB 
+Jabs (1 + cos ip) /2 (o+ k a 2k + h^a^ - &2k a 2-k - & 2ka2-k) 7kABi 

using the lattice structure factors ^aa — cos(k x /2) cos(fc y /2), j^bb — {cosk x + 
cosk y )/2 and J^ab = cos(k x ^/2). 

As stated, we use the coupled cluster method (CCM) to find the exact 
ground state of (||). To do this we notice the following property of H 

£if k = £(/f k + #_ k )/2 = £# k ; =*[ff£,i£,]_=0 Vk,k'. (4) 

k k k 

Hence it is possible to treat each seperately within the CCM, since they all 
commute with each other. So we have to deal with a bosonic system with eight 
different bosonic operators ai± k , a2±k, &i±k, &2±k denoted with a%, . . . , as- 

The ket and bra ground state of such a system (i.e a many-mode bosonic 
field theory with bosonic operators a,;, a+ in the Hamiltonian) in CCM-SUBZ 
approximation is given by [|5|, 

W=e s |0), S = Ei 1 ,i 2 ,..., il ^ 1 , i2 ,..,i I «---< ) ^ 
= <0|5e- s , S=l + }2 h - 2 A il , i2 ,..., i! a il a i2 • • • a 4i ' 1 ) 
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where |0) is the bosonic vacuum state (i.e. aj|0) = 0), and A il ... and A il ... are the 
CCM correlation coefficients. These coefficients are calculated by two systems of 
equations (one of them is a system of nonlinear equations). 

f) ff f)FT 

' =0, j- A -0, H={*\H\9), (6) 



dAi^ ' dA lt ... H 

using the expectation value (H) of the Hamiltonian, i.e. the ground state energy. 

Note, that the CCM-SUB2 approximation (i.e. having only quadratic terms 
of bosonic operators in S and S (^J)) gives the exact ground state of a quadratic 
bosonic Hamiltonian, since the ground state wave function of such a Hamiltonian 
has the form |^) = expfj^- fijafa~j~]\0), which can easily be shown using a Bo- 
goliubov transformation (see appendix). Therefore the CCM correlation operator 
S (and S respectively) (|^) consist of products of two bosonic creation operators 
only, all other coefficients Ai ll ... i l are zero; so we just have to use SUB2. 

To calculate the CCM equations (ra) easily using computer algebra, we make 
use of the Bargmann representation [Q, f| 



&z, a |0> <=> 1, (0\f(a,a+)\0)& f(-f,z) 

az dz 



(7) 



which maps the original many-mode bosonic field theory into the corresponding 
(classical) field theory of complex functions in a particular normed space. So 
instead of bosonic operators we just have to handle with (complex) numbers and 
differential operators, which is much easier. Once the (partial nonlinear) equations 
are obtained they can be solved numerically. 

3. Results and conclusions 

We apply the CCM-scheme described above to calculate the exact ground 
state of (||) and by doing this getting a spin wave approximate ground state of 
the model ([!]). We discuss the energy as a function of the angle between spins 
of the two subsystems A and B and the correlation between spins of different 
subsystems as a function of Jab (Fig-0)- We find as a typical order from disorder 
effect, that the degeneracy of the ground state with respect to the angle 99 is lifted 
by quantum fluctuations and a collinear ordering (ip — 0, it) is stabilized. This 
can clearly be seen by the energy vs. ip picture in Fig.|l| and by the correlation 
(SiSj)A,B vs. Jab, which is zero in the classical case, independent of the value of 
Jab (for \Jab\ < %V JaaJbb)- In the quantum case however that correlation does 
depend on Jab, showing again an order effect induced by quantum fluctuations. 

In addition we find a lowering of the magnetic order within the subsystems 
A and particular B by frustrating Jab in the quantum case. 
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A. Proof that CCM-SUB2 gives exact ground state 

Using the fact, that a Bogoliubov transformation (3 V = X] At ( M ^i/ a M ~ v *nv a t,) 
exactly diagonalize a quadratic bosonic Hamiltonian with the bosonic operators 
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Fig. 1. Ground state of the model (|l|): (a) energy as a function of the (classical) angle tp 
(in units of n) between the two subsystems A and B (note, that in the classical case the 
energy does not depend on cp); and (b) correlation between spins of these two subsystems 
in dependence on Jab- 

di, af, one can easily show that its ground state must have the form |^) = 
|0), by showing that = W. We use the Bargmann 



,+„+ 



ex p[Eij fi3»i « 

representation 



(0) and 



get 



exp 



E^ fij z i z j 



= Vz 7 - 



fj. i fi 

and this last matrix equation is allways fulfilled for some fi t 
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